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Abstract 

It is demonstrated in the context of the simple one-dimensional example 
of a barrier in an infinite well, that highly complex behavior of the time 
evolution of a wave function is associated with the almost degeneracy of levels 
in the process of tunneling. Degenerate conditions are obtained by shifting 
the position of the barrier. The complexity strength depends on the number 
of almost degenerate levels which depend on geometrical symmetry. The 
presence of complex behavior is studied to establish correlation with spectral 
degeneracy. 
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Tunneling processes have become of considerable interest as one of the possible mecha- 
nisms for creating highly complex behavior in the structure of the quantum wave function. 
Tomsovic and Ullmo ||T| found that there is an interesting correlation between classical 
chaotic behavior and the rate of tunneling in the corresponding quantum system. The 
conclusion of their study is that chaos facilitates tunneling. 

On the other hand, Pattanayak and Schieve found chaotic behavior in the semi- 
classical phase space (defined by expectation value) of a one-dimensional time-independent 
Duffing oscillator where new variables, associated with dispersion of the quantum states, are 
defined and included in the description of the system. They concluded that quantum tun- 
neling plays a crucial role for the chaotic behavior in the corresponding semi-classical maps. 
They have argued that the spectrum becomes more complicated in the neighborhood of 
the separatrix. 

In a recent study, we have considered a model in which tunneling leads to highly com- 
plex behavior of the quantum wave function and its time dependence. The spectrum, as 
anticipated by Pattanayak and Schieve [Q, indeed makes a transition to more complex 
behavior in the presence of a classical separatrix [Q. It is clear that the almost degeneracy 
of levels is necessary for the existence of significant tunneling. We directly investigate, in 
this work, the effect of almost degeneracy in the presence of tunneling on the complexity 
and the behavior of the development of the wave function. This criterion is, in fact, closely 
analogous to the criterion of overlapping resonances for the onset of classical chaos 0. 

The model we shall use is related to the one we previously explored P], i.e., a barrier 
embedded in an infinite well. By displacing the barrier in the double well system (to the 
right or to the left), certain positions are passed where the system becomes strongly almost- 
degenerate. These positions occur at almost commensurate intervals. It is exactly for those 
positions that one may find significant tunneling accompanied by complex behavior. We 
show, moreover, that in the cases of very high degeneracy, tunneling from left to right 
has exponential decay, on a significant interval of time, but at other positions, where almost 
degenerate conditions are somewhat weaker, the transition curve develops strong oscillations. 
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In a study by Nieto et al. |0, it was shown that tunnehng in an asymmetric double well 
is a very sensitive function of the potential. The behavior of the development of the wave 
function under these conditions was not, however, discussed there. 

The calculation in this work is done for a square barrier of height V = b and width 
w = X2 — Xi = 2 (xi is the left boundary of the barrier and X2 is the right boundary; 
we take % = 1,2m = 1) embedded in an infinite well of width 21 = 110 (interval (—/,/)). 
In the calculation, an exact analytic expression is evaluated on the computer; there is no 
accumulation of error for large times, since there is no integration over time. 

We first discuss the energy spectrum according to the location of the barrier. In Fig. [l| 
one can see the lower energy levels "almost crossing" (the levels do not cross, but can become 
very close to each other) at several locations of the barrier. In the middle (the position c of 
the center of the barrier is taken at c = 0; generally, c = |(xi + X2)) every energy level is 
almost degenerate. In other discrete locations, one finds almost degeneracy for every second 
level, every third level, and so on. 

As an approximation to our model, consider two separate infinite wells with widths l + Xi 
and I — X2- The energy levels of the two separate wells exhibit very similar behavior to that 
of the finite barrier, but in this case, exact degeneracy occurs, according to the geometrical 
configuration of the system. For X2 = —xi, we have complete symmetry and all levels 
are degenerate. When the width of the left well is twice the width of the right well or 
vice versa, one obtains degeneracy for every second energy level (of three levels, two are 
degenerate). If the width of the left well is one third of that of the right well, every third 
level is degenerate, and so on. This follows from the relations En^ = h?{'KniY /2m{l + xiY 
and En^ = h^^-Krij-Y /2m{l — X2Y', if ni/rij. = (/ + Xi)/{1 — X2), then E^n = {ni and Ur 
are positive integers). The locations for which these degeneracies occur are : c = [{rir — 
ni) I [rir + ni)]{l — w/2). In every rii + rij. levels we have at least one degeneracy. In our model 
(barrier in infinite well), the behavior is very similar to the problem of separate wells for the 
lower energy levels (for the higher energy levels there are, in fact, no crossings); degeneracy 
locations are, moreover, shifted slightly forwards the center. 
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In order to study the tunneling process, we constructed a wave packet approximating the 
form ipi^x, 0) = cexp(— (x — xqY/ (4cr^) + ik.j.) (where a = 5, fc^ ~ 0.45 and xq = {xi — l)/2), 
from 28 or less of the first energy levels, located in the middle of the left side of the barrier 
(the normalization of the wave packet is approximately 0.9999). The average energy is 
approximately 0.1 x V. We then measured the maximum probability to be on right side of 
the barrier during a very long time interval {tmax ~ 2 x 5 x 10^). The calculation is done for 
the central region of the well (from — //5 to //5) to avoid a phase space imbalancing effect (if 
the barrier is located too close to the left side, for example, the probability to be in the left 
side is much smaller than the probability to be in the right side, just because the available 
space is much less). It can be seen clearly from Fig. that just several positions (which we 
call RE, for resonance enhancement) allow tunneling, while in most regions the wave packet 
is trapped in the left side. The strongest RE is found in the center of the well where we 
have complete symmetry. Fig. |^ and Fig. ^ show complete correspondence; the strength 
of the RE's depend on the number of degenerate energy levels (Fig.|^, inset). The second 
strongest RE in the picture is found where the system has almost degeneracy for every fifth 
level (c ~ —54/5 and c ~ 54/5). 

An additional factor that one must consider is the projection of the initial state (Gaussian 
wave packet) on different eigenfunctions, i.e., the coefficients of the representation. It is clear 
that the strength of an RE depends on the number of almost degenerate states that have a 
large overlap with the initial state. Thus, in order to get strong RE, there must exist at least 
one pair of eigenstates, |j — 1 >, |j >, which fulfill two conditions : 1) Almost degeneracy 
of levels {Ej_i ~ Ej), and, 2) the scalar product of the eigenstate with the initial state is 
large enough (i.e., in our case, {ip{0), appreciable compared to unity). Condition 1 forms 
a general underlying symmetry of the system, while condition 2 is a requirement for the 
symmetry effects to be realized. 

In fact, almost degeneracy of levels (condition 1), Ej^i ~ Ej, implies symmetry proper- 
ties of a pair of eigenstates, | j — 1 > and \j >. Let us denote the part of the eigenfunction 
on the left side of the barrier as 0^, and on the right side of the barrier as (pn (for simplicity 
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we neglect the function under the barrier, since the eigenf unctions are small in this region). 
Almost degeneracy of levels and orthogonality implies that the eigenfunctions are almost 
the same on one side the barrier, and opposite in sign on the other side of the barrier, i.e., 
(f^Lj-i ~ 0L,i = 4>L and 4>R,j-i ~ ~4>R,j = or vice versa. Moreover, the eigenvalue 

condition requires that : ~ /r while the normalization condition requires : 

Xl + Xr ~ 1- Thus, |0Lp ~ Ir\4>r? ~ I- These properties implies symmetry 
and antisymmetry in the central position. We conjecture that this symmetry is the essential 
property of the central position (c = 0). 

The wave function for the main, central RE exhibits a complex behavior for the evolution. 
This complexity is due to the large number of almost degenerate levels. When we measure 
a physical quantity, the difference between levels determines the time dependence (the time 
dependent phase is computed according to AEij). The very small frequency due to almost 
degeneracy implies a very large recurrence time. On the other hand, large energy differences 
are associated with short time scales. The influence of these types of frequencies can be seen 
in most of the results. However, this behavior does not occur for the total probability in the 
left side of the well as a function of time, as seen in Fig. |^a. The curves are smooth and 
do not reflect the influence of the short time scale. As explained before, at RE locations 
there exist, at least, one pair of eigenfunctions (pj-i, (pj which are approximately the same 
on the left side and f^} 4>j-i4>j is appreciable (approximately |). The influence of other 
(non- neighboring) eigenfunctions on Pie ft, tend to be small. The result is a combination of 
periodic functions (the number of these functions corresponds to the number of pairs that 
fulfil the two conditions for RE), with very small frequency differences. 

The transition from one side to another, when c = Cq = 0, exhibits approximate expo- 
nential behavior for times not too short or too long (Fig. ^). We observe similar behavior, 
but less clear, in other locations that have very strong RE (locations such as, c ~ —54/3 
and c ~ —54/2). It appears that this exponential decay is due to the behavior of a sum of 
periodic functions with very small different frequencies, as can be seen in Fig. |^a. After this 
interval of decay, Pieft enters a domain of large oscillations. 
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In Fig. we show also the results of the same calculation for some other RE locations. 
As expected, we find a periodic (or almost periodic) behavior. The Fourier transform (Fig. 
^) shows very strong frequency peaks that fit to the most dominant almost degenerate 
energy levels, and show clearly that very small frequencies dominate the motion. 

The almost degeneracy of levels that produces a high level of complexity i.e., chaotic-like 
behavior, occurs in the presence of strong tunneling. We have observed this effect in ref. 
0. We study here one of the most clear ways to display this connection. We compare, in 
Fig. I^a, results of four positions. For c\ = —10.89, there is a large RE with three pairs of 
dominant almost degenerate eigenstates, and for C2 = —3.63 there is a smaller RE with one 
pair of dominant almost degenerate eigenstates. The choice of C3 = — 1 results in no RE, no 
degeneracy and no tunneling at all. The main RE, at c = Cq = 0, reflects a very complex 
behavior as we have shown in ref. 0. 

In Fig. ^ the entropy deflned by S{t) = — J \ilj{x,t)\'^\'n\ilj{x,t)\'^dx is computed. For 
Co, shown in Fig. the entropy rises sharply accompanied by high frequency oscillations 
and then remains a long time in a "quasi-equilibrium state". The second RE, ci, shown in 
Fig. shows a tendency to recurrence after 280000 time units, while C2, shown in Fig. 
returns to almost the initial condition after approximately 80000 time units. The entropy 
for non-RE locations shows almost periodic behavior (Fig. the inset shows the structure 
at increased scale). A similar behavior can be seen for p{t) =<< > — < x >^>. 

Comparison between different locations shows that one can characterize the behavior by 
two time scales. The flrst is the time of approach to equilibrium (Atgg), and the second 
is the recurrence time (At^c)- The approach to equilibrium time corresponds to averaging 
small frequencies (i.e., < AE > over all energies that satisfy almost degeneracy). One can 
easily identify At^q from Fig. ^a (for cq) and from Fig. while Atrc can be calculated 
analytically from the known eigenvalues. It appears that can give a measure for the 
complexity of the behavior of the system, as seen in these results. In all cases the, Atgq, is 
approximately the same, while Atrc is changed drastically from cq to C3 (one can not recognize 
recurrence in the Cq location, while Ci and C2 exhibit almost recurrence as mentioned in 
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previous paragraph). Thus, the ratio of these two time scales is largest for C3 (largest RE 
and maximum complexity), and decreases when the RE's become stronger (or when the 
complexity becomes stronger). 

In Fig. ^, we compare the < p >.< x > maps, sampled at peak times (times at which 
the wave function forms peaks). For each location of the barrier (co,Ci,...), we measure 
the time between peaks of \ip{x 1,1)1"^. At cq (a) the map moves toward the center and then 
accumulates in the central region. The map of ci (b) shows some ordered lines, and it is easy 
to see that the wave packet stays most of the time in the left side. More ordered behavior 
appears in the map of C2 (c), where the lines appear as semi-periodic paths. [The mapping 
oscillates back and forth, each time adding points to paths.] Finally, C3 (d) provides ordered 
maps, as expected. The same pattern can be seen also in < tt >:< p > (tt = dp/dt), at 
peak times, and from the Poincare maps < p >.< x > sampled when p{t) is minimum (i.e. 
71 (t) =0 and dn/dt > 0) [0. A period doubling behavior (small circles within large circles) 
can be seen in the < p >:< x > plane, and in the < tt >:< p > plane, as shown in ref. 0. 
The behavior becomes more ordered, as the RE height decreases, while for C3 there are just 
large circles. 

We have shown in this study that tunneling in the presence of almost degeneracy of levels 
provides necessary and sufficient conditions for the development of complex behavior of the 
wave functions of a quantum system. A large number of almost degenerate levels induces a 
high level of complexity. In general a pair of almost degenerate levels has a pair of equivalent 
eigenfunctions for which one is non-alternating and the other alternating. Moreover, in this 
case, the probability to be in the left side is approximately equal to the probability to be 
in the right side, without connection to the barrier position. These conditions suggest an 
underlying symmetry property which depends strongly on the geometry of the system. We 
wish to thank E. Eisenberg, M. Lewkowicz, I. Dana, and R. Berkovits for many useful 
discussions. 
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FIGURES 



FIG. 1. Square root of energy eigenvalues as a function of the center of the barrier. 

FIG. 2. Maximum probabihty to be in right side of the barrier. The inset shows partial 
correspondence between RE (resonance enhancement) and almost degenerate levels. 

FIG. 3. (a) Probability to be in left side of the barrier as a function of time for three largest 
RE's (in the interval (—1/5,1/5)). (b) The power spectrum of (a). 

FIG. 4. The entropy function as a function of time for different locations, (a) c = cq = (b) 
c = ci = —10.89 (c) c = C2 = —3.63 (d) c = C3 = —1, the inset shows an enlargement of typical 
periodic oscillations. 

FIG. 5. The accumulation of points in the < x >:< p > plane according to peak times for 
different locations, (a) c = cq = (b) c = ci = —10.89 (c) c = C2 = —3.63 (d) c = C3 = —1. The 
doted line indicates the position of the center of the barrier. 
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